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Abstract 

Quantum Energy Inequalities (QEIs) are results which limit the extent to which 
the smeared renormalised energy density of the quantum field can be negative, 
when averaged along a timelike curve or over a more general timelike submanifold 
in spacetime. On globally hyperbolic spacetimes the minimally-coupled massive 
quantum Klein-Gordon field is known to obey a 'difference' QEI that depends on 
a reference state chosen arbitrarily from the class of Hadamard states. In many 
spacetimes of interest this bound cannot be evaluated explicitly. In this paper we 
obtain the first 'absolute' QEI for the minimally-coupled massive quantum Klein- 
Gordon field on four dimensional globally hyperbolic spacetimes; that is, a bound 
which depends only on the local geometry. The argument is an adaptation of that 
used to prove the difference QEI and utilises the Sobolev wave-front set to give a 
complete characterisation of the singularities of the Hadamard series. Moreover, 
the bound is explicit and can be formulated covariantly under additional (general) 
conditions. We also generalise our results to incorporate adiabatic states. 

Dedicated to Klaus Fredenhagen on the occasion of his 60th birthday. 



1 Introduction 



The classical minimally coupled scalar field, like most matter models studied in classical 
general relativity, obeys the weak energy condition (WEC). That is, the stress energy 
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tensor Tab obeys the inequality TabV°'v^ > for all timelike vector fields v"', which entails 
that observers encounter only non-negative energy densities. However, it has been known 
since 1965 that no Wightman quantum field theory can obey the weak energy condition 
[H], (see [ini [IZj for simple arguments as to why this is true). Moreover, under many 
circumstances there is no lower bound to the energy densities available in quantum 
field theory (QFT). This surprising feature of QFT has often been used to support 
proposals for exotic spacetimes, such as warp drive and traversable worm holes, which 
require WEC-violating matter distributions. In addition, the validity of the second law 
of thermodynamics is called into question by WEC violations in QFT [26]. 

With these concerns in mind, substantial effort has been directed to understand the 
magnitude and extent of negative energy densities permitted by QFT, starting with the 
work of Ford in 1978 [26]. Given the failure of pointwise energy inequalities, attention 
has been focussed on averages of the stress tensor along timelike worldlines or over 
spacetime regions. In many QFT models such averages turn out to obey Quantum 
Energy Inequalities (QEIs); that is, their expectation values are bounded from below 
as the state varies within the class of physically reasonable states [HI ttUl [HI HH [2Sl EZl 
fIE\ m] . Since their inception QEIs have been applied to a variety of physical problems; 
they form the basis of the arguments constraining exotic spacetimes such as the warp 
drive [21 EH] and traversable wormholes [271 EH] • 

In their most common form, QEIs place lower bounds on the expectation value of 
the averaged stress energy tensor relative to that obtained in a reference state. For this 
reason they are called difference QEIs. To give a specific example [TU], consider the 
minimally coupled scalar field of mass > in a globally hyperbolic spacetime {Ai,g) 
of dimension m > 2, and let 7 : M — be a smooth timelike curve (not necessarily a 
geodesic) with velocity v"'. Then for any real- valued / G C^(M) the QEI 

dtfit) ((.;VTr).(7(t)) - {v^vXrum)^ > -bd (i) 

holds for all Hadamard states uj, where the bound 
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(2) 



depends only on /, 7, and the reference state ojq (which may be any Hadamard state); 
note that it does not depend on the state of interest oj. Here the hat denotes the 
Fourier transform given, in our conventions, by f{C) = Jj^nd^x f{x)e^^'^. The quan- 
tity {v°'v^'T^I}^^)^q{x,x') is the (unrenormalised) point split energy density defined, in a 
neighbourhood of 7, by 

K^''T:r)c.o(^,^') = ( ^ yyiVa ® e^v,. + ^/.^i ® i ]a.,{x,x') (3) 



1 1 \ 

a=0 ^ 



where {e^}a=o,i,2,3 is a tetrad field satisfying eg|^ = (see Section 3 of [lOj for a 
more detailed discussion) and A^^^ is the two point function of the state cuq. Finally, 
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^9*(f'^f*'T^'f*)^Q denotes the pull-back 

^*{v^v''T:fX{r, r') = {v^v''T:fXMr)Mr')) (4) 

which may be defined rigorously as a distribution on using the techniques of microlo- 
cal analysis, which also guarantee that the bound (I2l) is finite. Similar bounds also hold 
for the free spin- 1/2 and spin-1 field in comparable generality and rigour [3 [121 IT3l fli] . 
Recently, quantum energy inequalities have also been proven for free spin-3/2 fields in 
Minkowski spacetime [IHl EZ] • 

As already mentioned, one application of the QEIs is to place constraints on ex- 
otic spacetimes. However, the bound given above, while valid in any globally hyper- 
bolic curved spacetime, depends crucially on the choice of a reference state. Although 
Hadamard states exist on any globally hyperbolic spacetimes [291130]. closed form expres- 
sions for two point functions are known only in very special circumstances. For instance, 
no such expression is available for any Hadamard state on the warp drive spacetime. Typ- 
ically these problems have been avoided by heuristic appeals to the equivalence principle 
to justify the use of Minkowski spacetime QEIs on sufficiently small scales. To date, this 
approach, while physically reasonable, lacks full mathematical justification and control 
over the scales on which it is valid. It would clearly be preferable to employ a lower 
bound which did not require the specification of a reference state and placed constraints 
directly on {T^^)^. Bounds of this type, known as absolute QEIs, have been estabhshed 
in fiat spacetimes [21] but the only curved spacetime absolute QEI is that of Flanagan 
[25] (see also [IS, H?]) which applies to massless free fields in two-dimensional globally 
hyperbolic spacetimes. This approach relies on the conformal invariance of the theory 
and does not generalise to higher dimensions or non-zero mass. However, it does provide 
the basis for a QEI on arbitrary positive energy conformal field theories in Minkowski 
spacetime, including interacting examples [18] . 

In this paper we present the first absolute QEI applicable to the scalar field of mass 
/X > in four- dimensions, by refining and modifying the argument presented in [10] . For 
averaging along timelike worldline 7, our result takes the form 

<ltf{t){v'^v'TX).{t)>-BA, (5) 



where 



[—(,,0 + "local curvature terms" (6) 



and T^^^^^H is constructed in the same fashion as {v'^v^'T^}^^)^^ but with (essentially) the 
first few terms of the Hadamard series replacing the reference two-point function. At 
the technical level, we invoke a refined version of microlocal analysis which keeps track 
of the order of singularities. 

The structure of this paper is as follows. In Section[2]we review the algebraic formula- 
tion of quantum field theory in curved spacetimes, review two (equivalent) formulations 
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of the Hadamard condition and give a detailed analysis of the singularity structure of 
the Hadamard series in terms of Sobolev wave-front sets. Section [3] contains our main 
result, theorem 13. which is then used to give a number of examples of absolute QEIs. 
Although our bound depends on a choice of coordinates, we describe how the depen- 
dence can be eliminated by restricting the choice of smearing tensor, thus providing a 
covariant formulation of our bounds. 

2 Quantum field theory in curved spacetime 

2.1 The algebra of observables and the first definition of 
Hadamard states 

We shall employ the algebraic framework for describing the scalar quantum field in a clas- 
sical curved four-dimensional spacetime {A4,g). Here is a four-dimensional smooth 
manifold (assumed Hausdorff, paracompact and without boundary) with a Lorentz met- 
ric gab of signature (H ). Furthermore, we require {J\4,g) to be globally hyperbolic, 

that is Ai contains a Cauchy surface. Where index notation is used, Latin indices will 
run over the range 0,1,2,3 unless explicitly stated otherwise, while Greek characters 
will denote frame indices and also run over 0,1,2,3 unless explicitly stated otherwise. 
We employ units in which c = h = 1. 

The minimally coupled scalar field (p obeys the Klein-Gordon equation (V^+/i^)0 = 0, 
where = g"'^'VaVb and > is the mass of the field quanta. Global hyperbolicity 
entails the existence of unique global advanced {E~) and retarded (E^) Green functions 

: C^{M) C°°{M) for the Klein-Gordon equation obeying 

{V' + ix')E^f = E^iV' + fx')f = f, (7) 

and 

supp E^f C J±(supp /) (8) 

for all / G C^{A4), where ^^(5*) denote the causal future (+) and past (— ) of a set 
S. One may use the set of smooth functions having compact support in Ai, C^{M.), 
to label a set of abstract objects {4>{f) \ f G C^{A4)} which generate a free unital 
*-algebra 2t over C. The algebra of smeared fields 'Ql{Ai,g) is defined to be the quotient 
of 21 by the following relations: 

i) Hermiticity, = 0(7) V/ G C^{M); 

ii) Linearity, 0(a/ + Pf) = a0(/) + /30(/') Va, /3 G C and V/, /' G C^iM); 

iii) Field equation, ^((V^ + /i^)/) = V/ G Co°°(A^); 

iv) Canonical commutation relations, [4>{f), </>(/')] = iE{f, f')l V/, /' G C^{M). 
Here, E = E~ — E^ is the advanced-minus-retarded Green's function for the Klein- 
Gordon operator and by E{f, /') we mean 

E{fJ')= [ dYo\{x) f{x){Ef){x). (9) 
Jm 
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It is relation (iv) that quantises the field theory. 

In this framework, a state is a linear functional u on 2t(A^, g) which is normalised so 
that a;(l) = 1 and is positive in the sense that lj{A*A) > for all A G 2l(A^, The two 
point function associated with the state a; is a bilinear map A^^ : C^{A4)®C^{J^) C 
given by A^(/, /') = Lj{(f){f)(f){f')). We will only consider states for which A^ is a 
distribution, i.e., A^^ G V{Ai x Ai). It is clear from (iv) that the antisymmetric part of 

1 Z' . . „ . , „, „.\ i 



2^A^(/,/')-A^(/',/)J =-E(/,/'), (10) 
is state independent. 

As already mentioned, we will largely be concerned with Hadamard states. There are 
two equivalent formulations of the Hadamard condition, both of which will be used in 
the sequel. The original definition, given in a precise form by Kay & Wald [32] , involves 
a local series expansion of the two point function A^ associated with a state u and is 
based upon Hadmard's work on the fundamental solution for hyperbolic operators. 

In order to give the precise formulation of the Hadamard series construction we first 
introduce some geometrical structures, following [39], |l2l HI] . We denote by X C x 
the set 

X = {{x,x') E Ai X Ai \x,x' causally related and 

J^{x) n J^{x') & J^(x) n J^{x') are contained 

within a convex normal neighbourhood} . (11) 

For each (x, x') G X let U^^x' be any convex normal neighbourhood containing J^{x)r\ J~{x 
and J~(x) n J^{x'). Then (cf. Lemma 3.1 in [l2j|) X = U(x x')e3e ^2^.2^' ^ open 
neighbourhood of X in Al x 7W on which the signec0 squared geodesic separation of 
points a is well-defined and smooth, and on which the Hadamard construction (to be 
described shortly) can be carried out. Any open neighbourhood of X defined in this way 
will be called a regular domain. 

For each A; = 0, 1, 2, . . . , we may define a distribution Hk G V{X) by 

1 f A^x,x') A a^{x,x') ( a+{x,x') \ 

where we have introduced a length scale ^ to make a/ ^^ dimensionles^ and the coefficient 
functions A, Vj and Wj will be explained below. We also set = A^/^/(4'7r^cr+). By 



^We adopt the convention that a{x,x') > if x,x' are spacehke separated, a{x,x') < if x,x' are 
timehke separated and a{x,x') = if they are null separated. In Minkowski spacetime, (R'*,??), for 
example a{x,x') = —Tjabix — x')°'{x — x')^. 

different choice of length scale £' can be absorbed into a redefinition of the local curvature terms 
Cab appearing in the renormalisation of the stress-energy tensor; see Section [ 
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F{a+), for some function F, we mean the distributional limit 



= lim FK), (13) 

where ae{x,x') = a{x,x') + 2ie{t{x) — t{x')) + and t is a time function; that is, V"t 
is a normalised future directed timelike vector field on X. We shall occasionally use the 
notation t = t{x) and t' = t{x'). The function A G C°°{X) is the van Vleck-Morette 
determinant bi-scalar and is given by 

V-^(a;)v-^(a;') 

The functions Vj and are found by fixing x' and applying (V^ + /x^) (8> 1 to if^ and 
equating all the coefficients of l/o"+, l/o"+, lno"+ etc to zero. This determines a system 
of equations (known as the Hadamard recursion relations, given in appendix which 
can be solved uniquely (in X) for the Vj series. The Wj series is specified once the 
value of Wo is fixed; we adopt Wald's prescription that wq = [18]. We remark that 
the k ^ oo limit of the right-hand side of (fT2|) has a nonzero radius of convergence in 
analytic spacetimes, but not in general [32]. 

Let A/" be a causal normal neighbourhood of a Cauchy surface C [32]; that is, C is a 
Cauchy surface for M and every double-cone J~^{x) fl J^{y) with x,y E Af is contained 
in a convex normal neighbourhood of {A4,g) (see Lemma 2.2 in [39] for the existence 
of causal normal neighbourhoods). We may further choose an open neighbourhood X* 
of the set of pairs of causally related points in Af x Af whose closure is contained in 
X r\{J\f xAf) and a cut-off function x : A/" x A/" ^ [0, 1] so that 

X\x, = 1 and x\(AfxAf)\x = 0. (15) 

See Lemma 3.3 in [12] for the existence of X* and x with these properties. 

Given the above, a state u on ^{Jli,g) is said to be Hadamard if for each k E N 
there exists a F^. G C'^{N' x J\f) such that 

= xHk + Fk (16) 

in A/" X A/". We remark that this definition can be shown to be independent of the choices 
of C, Af, t, X, X and X, [391 m]. 

In the special case in which A^ is a convex normal neighbourhood, we note that Ai 
would be a causal normal neighbourhood of any of its Cauchy surfaces and we could take 
X* = X = A4 X A4 and x = 1, so (fT6i) becomes A^ = if^ -|- F^ and holds on the whole of 
Ai X Ai. In the general case, it is easy to see (e.g., using the microlocal characterisation 
of the Hadamard condition) that if uj is Hadamard then so is its restriction to any 
open globally hyperbolic subset of Ai, considered as a spacetime in its own right. Thus 
A^ — Hk is for all k on any set of the form U x U where f/ is a globally hyperbolic 
convex normal neighbourhood. As every point x E Ai has such a neighbourhood Ux we 
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may conclude that A^^ — Hk is for all k in an open neighbourhood of the diagonal of 
the form IJxeAi ^ ^^^^ refer to any such open neighbourhood as an ultra-regular 

domain. 

Note: The need to introduce the notion of an ultra-regular domain only came to light as 
the final version of this paper was prepared for publication, and after |45j had gone to 
press. One of us (CJS) would like to warn the reader that some results in hold on 
ultra-regular domains as opposed to the stated regular domain; with this modification, 
the results of ^5] are unchanged. 

2.2 Renormalisation of the stress tensor 

The Hadamard series construction forms the basis for the renormalisation of the stress- 
energy tensor in curved spacetimes, to which we now turn. The classical stress tensor 
of the real scalar field 

Tah{x) = (Va ® Vfe - ^gabg^^^c ® Vrf + ^yu'^abl ® 1^ ® ^) (x, x) (17) 

must be renormalised in QFT owing to the divergent behaviour of the two point function. 
Define the point-split stress-energy operator (which should not be confused with the 
stress-energy tensor itself) by 

^af* = V„ ® Vb' - Igab'g^^'Vc ® Vd' + ^fi^gab' 1 ® 1 (18) 

near the diagonal in A4 x A4, where gab'ix^x') is the parallel propagatoiH. If u is 
a Hadamard state we may define {T^b^)uj{x) at any point x G by the following 
procedure: 

a) note that A^; — Hk G C^(X) for k >2 and any ultra-regular domain X, so T^^}^^[A^ — 
Hk) is defined and continuous near the diagonal in 7V1 x Ai; 

b) define 

{T^nM = lim <7,^'(x,x')r,\f*(A^ - H,)ix,x') (19) 

for k>2; 

c) make finite corrections to {T^i)u} in order to obtain a conserved tensor {Tll^)^^{x) with 
the correct properties in Minkowksi space. 

Step (c) is needed because the tensor (T^J^)^^ is not covariantly conserved and cannot 
be considered as an appropriate stress-energy tensor (it could not be inserted on the right 
hand side of the Einstein equations, for example). However, it turns out that V°(T^^'^)^ 
is of the form VhQ where Q is a local quantity, determined up to a constant [IS] : 
subtracting Qgab by hand from {T^{^)u} we therefore obtain a conserved quantity. The 
undetermined constant in Q is fixed by the requirement that in Minkowski spacetime the 

■^Thus, if G Tx'M.^ gab' (a;, x')v^ is its parallel transport to T^M along the unique geodesic joining 
x' to which is well-defined sufficiently close to the diagonal in x M.. 
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vacuum expectation value vanishes. If we require that the difference (T^l^)^ — {T^h^)^;^ 
should be given by 

(TT)^ - {TTU = lim g,'' {x, x')T:f' (A. - A.J (x, x') (20) 

then any remaining finite renormalisation must take the form of a state-independent 
conserved local curvature term Cab that vanishes in Minkowski space, and the finite 
renormalised expection value of the quantum stress energy is given by 

(TD.W = (rfir).(x) - Q{x)gab{x) + Cab{x) . (21) 

We take the view that the tensor Cab is a necessary part of the specification of a given 
species of scalar field, alongside the mass and curvature coupling. Given sufficient ex- 
perimental accuracy Cab should, in principle, be measurable. 

The renormalisation prescription we have outlined above is vulnerable to the criticism 
that the Qgab term needed to restore conservation of {TH^)^ is only found to be a 
local curvature term a posteriori. Moretti [ID] has shown that this problem can be 
circumvented by an alternative construction of the quantum stress tensor. The basic 
idea is to modify the classical stress energy tensor Tab by adding a term of the form 
aipiy"^ + yU^)v5 for constant a. While this addition does not affect the classical physics, 
it has a non-trivial quantisation. A judicious choice of a ensures that the quantised 
stress energy tensor is conserved a priori (see theorem 2.1 of [10]) and agrees with the 
usual quantisation up to conserved local curvature terms. Although Moretti's approach 
is certainly elegant, it turns out that the usual quantisation is better adapted to the 
derivation of QEIs. In particular, our argument relies crucially on being able to write 
T^^/'* in a symmetric form which is not possible with a term of the form 1 ® (V^ -f jj?). A 
similar problem arises for the non-minimally coupled scalar field. In this case one must 
smear the stress-energy tensor even to obtain an inequality on the classical field [21], 
which necessitates a more complicated analysis at the quantum level [22] . 

2.3 The wave-front set and second definition of Hadamard state 

The above discussion shows that Hadamard states are characterised by their singularity 
structure. For this reason the techniques of microlocal analysis, which focus attention on 
singular behaviour, are ideally suited to this theory. This realisation has led to a number 
of important developments in the theory of quantum fields in curved backgrounds, follow- 
ing initial work of Radzikowski [12] , particularly in regard to renormalisation [3], |33l [3l] . 
In addition, the theory of (smooth) wave-front sets is a key tool in the proof of general 
difference QEIs [TOl [121 El E] in general globally hyperbolic spacetimes. Our absolute 
QEIs require the finer control on singularities of distributions afforded by the Sobolev 
wave-front set. In this subsection we briefly review the deflnition of the smooth and 
Sobolev wave-front sets and explain how they may be used to give a purely microlocal 
deflnition of the Hadamard condition, as flrst identifled by Radzikowksi [12]. In addi- 
tion, we will state a result of Junker and Schrohe [38] on the Sobolev wave-front set of 
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the two-point functions of Hadamard states. This will form the basis of our analysis 
of the Sobolev wave-front sets of individual terms in the Hadamard series in the next 
subsection. 

To begin, let u G T>'{W^) be any distribution. We say that u is smooth at x' if there 
exists an open neighbourhood O C M"^ of x' and a smooth function ip G C°°((9) such 
that u{f) = Jjg„ d"^x ip{x)f{x) for all test functions / G C^{0). The singular support, 
singsupp M, of a distribution u G V{W^) is the complement in R"* of the set of all points 
at which u is smooth. In particular, a distribution is smooth if and only if its singular 
support is empty. 

While the singular support tells us 'where' a distribution u fails to be smooth, Fourier 
transforms of localisations of u contain additional information. A covector ( G \ {0} 
is a direction of rapid decay for u at x ii there exists a conic neighbourhood T C M"* \ {0} 
of ( and a localiser x ^ C^(R™) which does not vanish at x such that 

(1 + 1^1)^1x^(01 — ^0 as^^ooinr, VA^ G N . (22) 

The set of singular directions of u E ©'(M"^) at x, S^(m), is the complement in \ {0} 
of the set of directions of rapid decay of u at x. The wave-front set of u assembles this 
information in a convenient way (see Section 8.1 of [35] for more detail). 

Definition. The (smooth) wave-front set WF{u) of a distribution u G X''(M'") is 

WF{u) = {(x,0 G M'" X (M- \ {0}) I i G S,.(m)} (23) 

As an example, it is easy to verify that the Dirac 5 and Heaviside 9 distributions 
have the following wave-front sets. 

WF{5) = WF{e) = {(0, I ^ e K \ {0}} . (24) 

The wave-front set is a closed cone in M"^ x (M™ \ {0}), whose elements transform as 
covectors under coordinate transformations (see Theorem 8.2.4 in [35]). Accordingly, 
the definition of wave-front set may be extended to distributions on a smooth manifold 
Ml in the following way: We say (x, ^ WF{u) C T*J^ \ {0} if and only if there exists 
a chart neighbourhood {k,U) of x such that the corresponding coordinate expression of 
(x,0 belongs to WF{u o k^^) C x (W" \ {0}), where m is the dimension of ^A. 

The Sobolev wave-front set provides greater structure on the information in the 
wave-front set. Recall that the Sobolev space H^{W^), s G M, is the set of all tempered 
distributions u on such that 

/ d™e(i+Kr)is(e)r<oo. (25) 

We summarise some relevant properties of Sobolev spaces for convenience. 
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Proposition 2.1. The Sobolev spaces H^{M."^) have the following properties: 



i) Ifs>k + m/2 then if ^(M™) C C"=(R") for keN; 

ii) H'iR"') C iJ"'(M")Vs > s'] 

Hi) if*(M™') is closed under multiplication by smooth functions. 



Associated with tlie scale of Sobolev spaces, there is a refined notion of the wave-front 
set. Just as WF{u) informs us where a distribution u fails to be smooth, the Sobolev 
wave-front set WF^{u) contains information about where in phase space the distribution 
fails to be H^. 

Definition. The distribution u G ViW^) is said to be microlocally at {x,S,) € 
X (M™ \ {0}) if there exists a conic neighbourhood T of ^ and a smooth function 
cp e C^{W), (p{x) ^ 0, such that 



The Sobolev wave-front set WF'^{u) of a distribution u G T>'{W^) is the complement, in 
T*W^ \ {0}, of the set of all pairs (x, ^) at which u is microlocally H'^. 

It is easy to verify, for example, that 



and we see how this refines the information in (1241) . Like the wave- front set, WF'^lu) 
is a closed cone in T*M" \ {0}. Furthermore, part (ii) of propostion 12.11 entails that 
WF'ju) C W F' (u) C WF{u) for all s < s'. In fact, one may show that WF{u) = 
UsgK ^-^*(^)- Additionally, if ^9 G C^(R"') does not vanish in a neighbourhood of x 
then (x, ^) G WF^i^u) if and only if (x, ^) G WF'^^ipu); we also have WF^i^u + w) C 
WF'^{u) U WF'^lw). One may show that WF^{u) can be characterised in a coordinate- 
independent way as a subset of the cotangent bundle which then permits the definition 
to be extended to distributions on a manifold by referring back to (any choice of) local 
coordinates (see, e.g., remark (i) following Prop. B.3 in [3H]). We shall occasionally use 
the notation u G H(^^{M.) if WF^{u) = for a distribution u G V'{M) (see also the 
remarks following definition 8.2.5 of [36]). 

In jl2] , Radzikowski proved the remarkable result that the definition of a Hadamard 
state in terms of the series construction previously given is equivalent to a condition on 
the wave-front set of the two-point function. Namely, the wave-front set is required to 
lie in a particular subset of the bicharacteristic set of the Klein-Gordon operator, which 
we now define. 



Denote by 7^ = {(x,0 e T*M \ g''\x)iaib = , ^ 0} the set of nonzero null 
covectors over Jvl. Since (A^, (?) is time orientable we may decompose TZ into two disjoint 
sets defined by = {(x,0 e 7^ | ±^ > 0} where by ^ > (^ G T*M) we mean 




(26) 



WF'^\e) = WF'{6) = 



{(0,0 kGM\{0}} s>-l/2 
s<-l/2 



(27) 
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that is in the dual of the future light cone at x. We define the notation (x, ^) ~ (x', 
to mean that there exists a null geodesic 7 : [0, 1] — > such that 7(0) = x, 7(1) = x' 
and = 7^(0)5'a6(x), = 'y^{l)gab{x'). In the instance where x = x', (x,^) ~ (x, 
shall mean that ^ = ^' is null. Then, the set 

C = {{x, e; x', e') G 7^ X 7^ I (x, ~ (x', 0} (28) 

is the bicharacteristic relation for the Klein-Gordon operator. We also define the related 
sets 

C+- = {(x,e;x',-neC|e>0} (29) 

and 

C-+ = {(x,-e;x',r)eC'|e>0}. (30) 

We may now state the relevant portion of Radzikowski's equivalence theorem [32 j : 

Theorem 2.2. Let {M.,g) be a four- dimensional globally hyperbolic spacetime and sup- 
pose A G V{A4 X A4) satisifies the Klein-Gordon equation and has antisymmetric 
part iE/2 modulo smooth functions. Choose a Cauchy hypersurface C, a causal normal 
neighbourhood N of C and a time function t. Then, the following two conditions are 
equivalent: 

i) A has the Hadamard series structure given by ( fl^) on Af x M , 

ii) WF{K) = C+-. 

As a consequence of this equivalence theorem we may now adopt the condition that 
WF{K^^) = C"^~ as the second definition of the state u, on ^{A4,g), being Hadamard. 

Junker and Schrohe [M] apphed the theory of Sobolev wave-front sets to study both 
Hadamard states and the larger class of adiabatic states (see §3.2.3p . In particular, 
lemma 5.2 of [38j gives the Sobolev singularity structure of the two point function of 
Hadamard states. 

Theorem 2.3. Let uj be a Hadamard state on %{M.,g) where M. is a smooth four- 
dimensional globally hyperbolic spacetime. Then, the two point function, h.^ eT)' { AA. x Ai), 
associated to uj has the following Sobolev wave-front set: 

WF'(K) = { l a . (31) 

2.4 Sobolev microlocal analysis of the Hadamard series 

We will now employ theorem 12.31 to study the Sobolev wave- front sets of the individual 
terms in the Hadamard series, working within an ultra-regular domain X on which the 
distributions l/cr+ and In a+ featuring in ( fT2l) may be defined. (The length scale i 
will be suppressed from now on.) We shall establish the following statement: 

WF^-'^^^anna^) C WF^{l/a^) = | ^ ^ . (32) 
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In order to do this we observe that the terms appearing in the Hadamard series are 
(loosely) related to one another via differentiation. If P is any partial differential oper- 
ator of order r on a smooth manifold Ai, i.e. in local coordinates 

P= J]pa(x)(-z9)" (33) 

|a|<r 

where a is a multi-index and pa are smooth functions, then the principal symbol, Pr{x,C,), 
of P is 

Prix,0= $^P«(x)r- (34) 

\a\=r 

The characteristic set of P, CharP, is the set of (x, ^) G T*Ai \ {0} at which the 
principal symbol vanishes. Corollaries 8.4.9-10 of encapsulate the effect of partial 
differential operators on the Sobolev wave-front set of a distribution: 

Lemma 2.4. Let M. he a smooth manifold. For u G V^M) and any partial differ- 
ential operator P of order r with smooth coefficients then WF'^{Pu) C WF'^^'''{u) and 
WP'+'^iu) C WF'{Pu) U CharP. 

Lemma 12.41 enables us to quantify our earlier observation about the relationship 
between l/o-+ and 0"-' lncr+: 

Proposition 2.5. Within an ultra-regular domain we have 

WF'+^+^{aHYia+) C WF'{l/a^) Vs G R andVj G {0} UN. (35) 

Proof. We employ induction on j. If i; G C'^{TM.) is a smooth vector field, then 
(f ■ V (S> 1) ln(T+ = [iy - V ® l)cr]/cr+. Hence, by lemma YTM, we have 

WF'+^{\Yi(y+) C WF'{l/a+) U Char {y-W®l). (36) 

As V is arbitrary, 

WF'+^{\na+) C WF'{l/a+) U [ p| Char (f ■ V ® 1) | (37) 

\v£C°°{TM) j 

and since 

Char {vV ®l) = {{x, x', i') G T*X \ | v{x) ■ ^ = 0} (38) 

it is clear that the intersection is empty and the statement holds for j = 0. Now suppose 
it holds for some j G {0} U N. The identity 

{v-V0 l)a^+^ \na+ = [{v V l)a] ((j + l)a^ \na+ + a^) (39) 

and the inductive hypothesis give 

WF'+^+^{a^+^ In a+) C WF'+^+^{a^ In a+) U Char {v V 1) (40) 
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and taking the intersection over all v G C°°{TAi) as before, we establish the result for 
j + 1 and hence all j G {0} U M by induction. □ 



We next prove the intuitively reasonable result that A^^ is as singular as the leading 
term in the Hadamard series. 

Proposition 2.6. Let u be a Hadamard state. Then, within any ultra-regular domain 
X, we have 

WF'{A^) = WF'{l/a+) VsgM. (41) 

Proof. Recall that for every A; G N there exists a G C^{X) such that A^^ = Hk + Fk- 
Hence, as WF'{anYia+) C W F'+i+^{anYia+) C WF'{l/a_^.), 

WF\K) C WF\A^ /a+)UWF\Fk) . (42) 

We remark that it is known (from, say, [42j) that A does not vanish where x, x' are null 
relatecEland as such WF^{K^) C W F''{1/ a+)\JW F''{Fk). Moreover, given any particular 
s we can always find a k sufficiently large such that WF'^{Fk) = and it remains to prove 
WF'{l/a+) C WF'{A^). Let {x,^;x',^') G WF'{l/a+) such that WF'-'{l/a+) = 
for any e > 0. Hence, by proposition 12.51 WF^{aHna+) C WF^~^~^{l/a+) = ^ for 
J > and we have - ^kj^-n'^a^ G H{^^{X). Therefore, {x,^;x',^') G WF%Hk) and 
by the nesting property (x, ^; x', G WF'^ (if^) for all s' > s. □ 

As a consequence of theorem 12.31 we now have the Sobolev wave-front sets of the 
constituent distributions in the Hadamard series. 

Corollary 2.7. The distributions 1/(T+, a-' ln(T_|_ G V{X), where X is an ultra-regular 
domain, have the following Sobolev wave-front sets: 

WF'{l/a+) = 

WF%anna+) C 

In consequence, we also have, for arbitrary j > —1, 

WF'iH,) c { %- I I (45) 

and 

WF-(H,,,-H,)c[%- (46) 

for j' > 0. 

^That is, A^/^ restricted to singsupp l/cr+ is non- vanishing. 





c+- 



s > 
s < 



1/2 
1/2 

s>j + 1/2 
s<J + 1/2 



(43) 
(44) 



13 



Proof. As we have already established, l/o"+ possesses the lowest order singularity which 
lemma 12.61 states is precisely that of A(^. The remaining results follow from proposi- 
tion ESI □ 



The remainder of this section is devoted to calculating the Sobolev wave-front sets 
of the advanced-minus-retarded fundamental solution E and a quantity if^ G V{X) {X 
an ultra-regular domain) defined by 

Hk{x, x') = i (^Hk{x, x') + Hk{x', x) + iE{x, x')^ , (47) 

which plays an important role in our main result theorem 13. II As iE is the antisymmetric 
part of A^, for all Hadamard u, theorem 12.31 implies that 

,yn>£)c{^""yc- ^ (48) 

Proposition 2.8. Within an ultra-regular domain, the Sobolev wave-front set of Hk 
satisfies 

' C^-\}C-^ s>k + ?,/2 
WF'{Hk)(l{ C+- s G [-1/2, A; + 3/2) . (49) 

s< -1/2 

Proof. Suppose first that s<fc + 3/2. It follows from the covariant commutation 
relations that, within an ultra-regular domain X, iE{x,x') = Hk^2{x,x') — Hi:^2{x',x) 
modulo C''+^(X). Hence there exists F G C''+^(X) such that 

[Hk - Hk]ix,x') = \Hk+2 - Hf,]{x,x') - - Hk]{x',x) + F{x,x'). (50) 

As A; + 2 > s, we have C^^'^{X) C H(^^{X), so all three terms on the right-hand side 
will belong to H{^^{X), using Eq. fHUl) as well. Thus Hk = Hk modulo H{^^{X) for all 
s < k -\- 3/2, which establishes (149 p for s in this range. For s > /c + 3/2 the result 
follows from the definition (jU]) of Hk, the wave- front set of Hk (and its behaviour under 
interchange of the arguments x and x') together with the rule for wave-front sets of sums 
of distributions. □ 

Finally, we end this discussion of the microlocal properties with the following result 
concerning the singularities of A^^ — Hk which follows directly from the proof of corollary 

Proposition 2.9. Within an ultra-regular domain, the Sobolev wave-front set of A^^ — Hk 

is given by 
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2.5 Restriction results and a point-splitting lemma 

In addition to the results of the previous subsection, our main result will make use of 
three additional technical results. The first, Beals' restriction theorem, enables us to 
restrict A^,Hk and their derivatives to certain submanifolds of x Ai. The second, 
taken from PU], shows that positive type is preserved under such restrictions, while the 
third result is a technical tool that enables us to write integrals over the diagonal on 
product manifolds in terms of their 'point-split' Fourier transforms. 

Our QEI results will encompass averages of the stress-energy tensor smeared over 
timelike submanifolds, e.g., timelike curves or hyperplanes, as well as averages over 
spacetime volumes. For these purposes, it is necessary to understand how restricting 
distributions (such as A^, and their derivatives) to a submanifold alters the Sobolev 
wave-front set. A theorem due to Beals (see Lemma 11.6.1 of [36j) tells us that, for 
suitably well behaved restrictions, the Sobolev order of the wave-front set is reduced 
by an amount proportional to the codimension of the restriction, while its elements are 
transformed according to the associated pull-back mapping. We will state a speciali- 
sation of Reals' result to the case we will need, in which we restrict from a product 
manifold 7W x to a submanifold E x E, where S is a submanifold of A4. Writing the 
embedding of E in as a map t : E — > and defining '(? = <.®<.:ExE-^A^xA1, 
the restriction of u G V{Ai x A1) to E x E may also be regarded as the formation of a 
pull-back ■})*u. Beals' result hinges on the relationship between the Sobolev wave-front 
set of u and the conormal bundle A^*E of E defined by 

A^*E = {(i(x), e T*M; a; e E, t*(0 = 0} . (52) 

Theorem 2.10 (Beals' restriction theorem). Let u G V\M. x M) and G C°°(E x 
E,7V1 X M) he defined as above, and suppose M. and E have dimensions m and n 
respectively. If {N*E x N*E) n WF%u) = for some s > m — n then '&*u is a well 
defined distribution in P'(E x E). Moreover, 

WF'-^"'-''\ru) C^*WF'{u) (53) 

where the set 'd*WF''{u) is defined to be 

i^*WF'{u) = {{t,L*{C);t',L*{0) e (T*E X T*E) I 

m,^;i{t'),a^WF%u)}. (54) 

The next result, theorem 2.2 of [TU], asserts that the positive type condition is pre- 
served under the restrictions carried out by Beals' theorem. 

Lemma 2.11. //, in addition to the hypotheses of Theorem \2.1(j\, u G X>'(A^ x M) is of 

positive type, then d*u is of positive type on E x E. 

Finally, we present a point-splitting identity for distributions of sufficient regularity. 
Beginning in x M", we have the following. 
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Lemma 2.12. For all u G Co(M" x R"), we have the identity 

[ d-tu{t,t)=lun [ ^e-l«l'^2(-e,0- (55) 

In particular, this holds if u e if'*(M" x M") fl x M") for s > n, by virtue of the 

Sobolev embedding of H'{W x M") zn ^(M" x M"). 

Proof. By definition of tlie Fourier transform, we fiave 



As tlie integrand is absolutely integrable on R^", Fubini's tlieorem permits us to reorder 
the integrations and perform the ^ integral first, thus obtaining 

= / dVdVV,(r-r')«(r,r') 



n (27r)" 



[ d''t(rt'ip,{t')u{t + t'/2,t-t'/2) 

d"tV,(t') / d''tu{t + t'/2,t-t'/2) (57) 



where (pe{t) = (47re)""/^e^'*'^''*^^^''. We have made the change of variables t = (r + t')/2, 
t' = T — t' (for which the Jacobian is unity) and reordered integrals using Fubini's the- 
orem again. As u is continuous and compactly supported, the inner integral exists for 
each t' and defines a continuous compactly supported function. The limit e 0+ exists 
and yields the value of this function at t' = because ipe is an approximate identity. 
This is the required result. □ 



Note that if ^ i— > u{—^,^) is absolutely integrable on M" then the dominated con- 
vergence theorem permits us to dispense with the limiting procedure on the right-hand 
sidell For our application, we will need a straightforward generalisation of the above to 
distributions on manifolds. 

Lemma 2.13. Let (S, /;,) be a n- dimensional pseudo-Riemannian manifold and u G 
£'{J1 X S) n H'^ijl X S) for s > n. Suppose the support of u is contained within U xU 
where U is a single coordinate chart of H with associated coordinate map k -.U —>■ M."'. 
Then 

[ dYol{x)u{x,x)= hm / ^e-^l«l'^(-e,0 (58) 

^Our hypotheses are strong enough to guarantee that u is absohitcly integrable on M" x M", and 
hence ^ ^ + + ^) is absolutely integrable a.e. in rj by Fubini's theorem. However, a simple 

proof of integrability for rj = was not forthcoming. 
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where U : n{U) x k{U) C is defined by 

U{x,x') = {\h^\^ (g) \h^\^u^){x,x') , (59) 
where = u o [k,~^ (g) k,~^) and is the determinant of the metric in coordinate chart 

K. 

Proof. We have 

/ dvol(x)-u(x,x) = / (rx\h^{x)\^^^u^{x,x) = / d"xf/(x,x). (60) 
Jt. Jr" Js." 

As hi^ is a positive smooth function bounded away from zero and therefore has smooth 
fractional powers, we may apply lemma 12.121 to the function U to obtain the desired 
result. □ 



3 An absolute quantum inequality 
3.1 Main result 

We now come to the statement and proof of our main result. Let S be any n-dimensional 
timelike submanifold of {A4,g) for 1 < n < 4, that is, h = L*g is a Lorentzian metric on 
E, where t : S — embeds the submanifold E in 7W. We also equip S with the time 
orientation induced from A4, so that non-zero future-directed causal covectors on (A^, (?) 
pull back to non-zero future-directed causal covectors on (T,,h). In our conventions a 
positive definite metric on a one- dimensional manifold is regarded as Lorentzian. As S is 
timelike, its tangent space TS can be annihilated only by covectors which can annihilate 
at least one nonzero timelike vector; in particular, all covectors in the conormal bundle 
N*I1 are spacelike. 

Our aim is to obtain lower bounds, as u varies among Hadamard states, on quantities 
of the form 

J dvol(x)/2(x)(g®g(A^-i/2))(x,x) (61) 

where Q = q°''W a + & is a partial differential operator with smooth real- valued coefficients 
and h defined on a neighbourhood of S and / G C^(S) is real valued. Note that 
A^ — H2 is so the coincidence limit is well defined. For simplicity it is convenient to 
assume in addition that S may be convered by a single coordinate chart with certain 
properties. 

Definition. A small sampling domain is an n-dimensional timelike submanifold S of 
{Ai,g) such that (ijH is contained in a globally hyperbolic convex normal neighbourhood 
in M.; (a) S may be covered by a single hyperbolic coordinate chart, i.e., a coordinate 
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system . . . ^ on S with d/dx^ future-pointing and timelike, and for which there 
exists a constant c > such that all causal covectors Ua on E obey 



c\uq\ > 



\ 



Y,n] (62) 

j=0 



(i.e. the coordinate speed of light is bounded from above). 

A sufficient condition for the existence of a maximum coordinate speed of light is that 
Hqq > e and | det(/ijj)"~i;^| > e for some e > 0. 

It is easy to verify (e.g. by using suitable normal coordinates) that every point of a 
general timelike submanifold S has a neighbourhood (in S) which is a small sampling 
domain. Thus any integral over a compact subset of a timelike submanifold may be 
decomposed into finitely many integrals over small sampling domains by a partition of 
unity. 

Suppose then, that S is a small sampling domain in {M..,g) with hyperbolic chart 
{a;"}a=o,...,n-i- We may express these coordinates by a map k : E — M*^, K{p) = 
{x^ {p) , . . . , x'"''^ {p)) and write = fi;(S). Any function F on S determines a func- 
tion = F o on S^; in particular, we have a smooth map : ^ Ai. The 
significance of k being hyperbolic is that the bundle TZ^ of (non-zero) future pointing 
null covectors on {A4,g) pulls back under so that 

t*7^+ C S« X r (63) 

where F C is the set of all Ua with uq > and satisfying ( l62i) . which means that F is 
a proper subset of the upper half-space x M"~^ of M" (here, we regard M"*" = [0, oo)). 

We now state our main result: 

Theorem 3.1. Let Tj be a small sampling domain of dimension n in {Ai,g) with hyper- 
bolic coordinate map k, and suppose Q is a partial differential operator of order at most 
one with smooth real-valued coefficents in a neighbourhood ofE. Set k = max{n-|- 3, 5}. 
For any real-valued f G C^(S) and any Hadamard state u we have the inequality 

J dvo\{x) f\x) {Q ® Q(A^ - H2)) {x, x)>-B> -00 (64) 



where 

d"e 



B = 2 



\h^\^fn®\h^\^f^r^{Q®QHk) 



(65) 



'iR+xii"-i (27r) 

is the determinant of the matrix n*h and t^iSxS— >-7Vlx7Vl is the map 'd{x, x') = 
{l (g) l){x, x'). 

Remarks: This bound depends nontrivially on the coordinates (and on any partition 
of unity used to reduce a general timelike submanifold into small sampling domains) . In 
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§3.31 we will discuss some classes of QEI averages which, in a sense, determine a natural 
choice of coordinates. Note that although H2 is sufficient to renormalise the left-hand 
side, the bound is given in terms of for k > b. Similar results hold when Q has 
higher order, for suitably modified values of k] we have restricted attention to the cases 
relevant to QEIs. 

Our strategy will be to mimic the proof of the general worldline quantum energy 
inequality presented in [10] but with Sobolev wave-front sets, as opposed to the smooth 
wave-front sets used in that paper. 

Proof of theorem \3.1[ We will break the proof into three parts. Part one will establish 
that 



dvol(x) f\x)Q ®Q{K- H2) (x, x) 



2 lim 



(27r)" 



\h^\-^f^®\K\-^fJlQ®Q{K-Hk) 



(66) 



for the given values of k. Part two contains a positivity result which enables us to discard 
the state dependent contribution to the right hand side of (l66l) to obtain the inequality 



dvol(x) f\x)Q ®Q{K~ H2) (x, x) 



> -2 lim / — ^ 

^^0+ Jk+xR"-i (2vr)" 



\K\'^f.®\K\-^f.r,{Q^QHj,) 



<,0 (67) 



Then, in part three, we show that the right-hand side of this expression is finite and 
equal to the required lower bound —B. 

Part one: We begin by observing that A^^ — H2 and A^^ — coincide on the diagonal 
in S X S, so we may write 



dvol(x) f\x) {Q ® g(A^ - H2)) (x, x) 

dvoi(x) f^x) r (Q ® g(A, - Hk)) (x, x). 



(68) 



The latter form has the merit that A^ — Hk is symmetric and more regular than A^ — H2. 
We have also written in the restriction map i}* explicitly, anticipating later steps in the 
proof. By hypothesis on k, we may choose s G (6, A; + 3/2), and Proposition 12.91 tells us 
that within an ultra-regular domain X G Ai x j\4 



As Q ® Q is at most second order, lemma 12.41 entails that 



(69) 



(70) 
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As the wave-front set WF'^ "^{Q ® Q{^lu — Hk)) is therefore empty and s — 2 > 4 — n, 
Beals' restriction theorem, theorem 12.101 yields 



so the point-sphtting identity, lemma 12.131 may be applied to give 



dvol(x) f\x) rq ®Q{A^- Hk) {x, x) 



(71) 



lim 



(27r)" 



(g) f^i^lQ ^ Q{A^ - Hk) 



(72) 



Then, as A^^ — Hk is symmetric and C'^, the integrand of fl72l) is invariant under ^ — ^, 



so we may replace the integration over R" with that over 



X 



pn— 1 



at the expense of 



a factor of two, thus obtaining ( !66l) . Note that we are now integrating over those ^ with 
^0 > 0. This particular half-space of R" is chosen because it contains the cone F defined 
after (1631). 



Part two: We now assert that the pull-backs of {Q (g) Q)A^ and {Q (g) Q)Hk exist 
separately, which enables us to split the integrand on the right-hand side of fl66|) into 
two parts. Moreover, we will show that the first of these, namely [ih^]^ ff^^lh^]^ f^'&l^Q'^ 
QA^'j^{—^,^), is nonnegative for all ^ G R". 

The existence of the required pull-backs follows because we have already observed 
that non-zero covectors in N*!! must be spacelike. As the covectors in the wave- front 
set of {Q®Q)A^ and {Q®Q)Hk are null, it follows that there is no intersection between 
their wave-front sets (at any Sobolev order) and N*Tj x A^*E . Therefore the pull-backs 
exist. 

To establish that [\h^\if^ \h^\if^§lQ ® QA^] (-^,0 > 0, we define a one pa- 
rameter family of functions /|(a;) = \hi^{x)\^ f i^{x)e^^'^ then, as G C^(R"), [I/^kI^/k ® 
® QAJ\ G £^'(Sk X Sk) where the Fourier transform is defined by m(C,^') = 
Vm G £'{W X R"). Therefore, 



\hn\~^f^®\h^\-^f^dlQ®QA, 



\K\~*f^®\K\-^f^dlQ®QA, 



<,0 



(ft /I) 



(73) 



(74) 



where we have exploited the fact that / is real valued. It is clear that if u is a distribution 
of positive type then Q ® Qu is also of positive type because Q has real coefficients. 
Accordingly, lemma [2.111 establishes that ■})*Q (8> QA^ is a distribution of positive type 
and the assertion is justified. Hence, we obtain the inequality ( 1671) provided that the 
limit on the right-hand side exists and is finite, which is the remaining step in the proof. 
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Note that the integral converges for each e > because the Fourier transform of any 
compactly supported distribution is polynomially bounded. 

A 



\K\-^f.(S)\K\-^f^i»lQ(S)QHk 



<,0 



Part three: Our aim is to show that J(^) := 

is absolutely integrable on the integration region M"*" x R"^^, for then we may conclude 
that the limit on the right-hand side of flUTj) exists by dominated convergence and equals 
—B (which is thereby finite). To do this, we introduce an arbitrary Hadamard state ujq 
and use the Hadamard series definition of Hadamard states to write = A^,, + Fj. for 
some Fk € C'^(X). We consider the contributions of these terms to I{C,) in turn. 

First, the results of Radzikowski and Beals entail that 

WF{^*,Q ® QA^J c KWFiQ ® QA^J c C^^" C C (^^ x , (75) 

where TZ"^ are the bundles of future- and past-directed null covectors defined earlier. 
Thus, we have 

WFi^lQ ® QA^J c 6:7^+ X tin- C (S. X F) X (S. x (-F)) (76) 

using equation ( l63l) and its obvious analogue for 71^ . 

By Prop. 8.1.3 in ^35j, it follows that the Fourier transform of localisations of ® 
QA^jy is of rapid decay outside the cone F x (— F); in particular we have rapid decay 
in the cone (— M"*" x M"^-*^) x (]R+ x R"^^) (here we have used the assumption that F 
is a proper subset of x M"^^ because k is hyperbolic). Accordingly we find that 
[I^kI^/k ® I/^kI^/k ® (5Aa,(,]^(— ^, is rapidly decaying in the integration region 
R"*" X ]R"~^ and is therefore absolutely integrable there. 

It remains to show that the F^ dependent contribution to /(^) is also absolutely 
integrable. As G C'^(X), we have {Q ® Q)Fk G C*^~^(X). Hence there is a constant 
c such that 



< (77) 

- (1 + + IfP)^'^"'^/' 



for all (^, ^') G X because is compactly supported. As (/c — 2) > n for the values 
of k given in the hypotheses, it follows in particular that left-hand side is absolutely 
integrable on M+ x M"^^. 

Accordingly, we have shown that / G L^(M+ x M"^^), and the dominated convergence 
argument mentioned above completes the proof. □ 



Two points should noted about the foregoing proof. First, the state ujq was introduced 
purely as a convenient way of showing that our bound is finite; the bound itself does 
not depend on any reference state. Second, in the difference QEIs studied in [TU] the 
Gaussian cut-off was not necessary, because the point-splitting lemma was applied to 
a smooth compactly supported function. Moreover, the place of Hk was taken by the 
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two-point function of a reference state A^^ and the fact that WF{Ai^g) = was 
used to show that the integrand decays rapidly in the integration region. This hne of 
argument was not available to us here, because (in contrast to Hk) has a portion of 
its wave-front set lying in (see Proposition 12.81) . 

In the next subsection, we will show how Theorem 13.11 may be used to obtain QEI 
bounds, by appropriate choices of the operator Q. 



3.2 Examples 

3.2.1 Worldvolume absolute quantum null energy inequality 

Our first example is a quantum null energy inequality (QNEI), that is, a lower bound 
on quantities of the form dvol {F°'^Tll^)^, where F""^ = n"'n^ and n°- is a smooth, 
compactly supported null vector field on {J^,g) that is future-directed where it is 
nonzero. We will show how Theorem 13.11 allows us to obtain an absolute QEI on 
Jj^ dvol {F"'^T^^)^. To do this, we suppose that F"* is supported within an open subset 
E that is a four-dimensional small sampling domain in {M.,g) with hyperbolic chart k. 
Noting that 



(F'^XD^(x) = hm (n'^Va ® n^'Vi^i " ^2) (x, x') + CabF''\ 



x), (7J 



we apply Theorem 13. II with Q = n-V and / G C^(E) chosen to be real- valued and to 
equal unity on the support of F"''. This yields the absolute QNEI: 



/ dvol (F'^XDc. > -2 / ^\K\^®\g.\^i{n-V®n.V)Hj 



A 

(-e,o 



+ / dvolF'^'^Cb (79) 
'm 



for all Hadamard states uo. 



Clearly the right-hand side of this inequality depends explicitly on the choice of 
coordinates k. In Section 13.31 we will explain how this problem may be removed by 
restricting the class of sampling tensors F""^ in such a way that there is a canonical class 
of coordinate systems, all of which give the same lower bound. 



3.2.2 Worldline absolute quantum weak energy inequality 

Our second example applies Theorem 13. II to the energy density sampled along a smooth 
timelike worldline 7. This was the situation studied in [lOj, where a difference QEI 
was obtained. We assume 7 is given in a proper time parameterisation as a smooth 
function 7 : J — M, where J is a possibly unbounded open interval of M and denote the 
four-velocity of the curve by m = 7. The curve forms a small sampling domain, with 
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the proper time parameterisation as a hyperbolic coordinate system, provided that the 
track of 7 can be contained in a globally hyperbolic convex normal neighbourhood in 
M. 

The classical energy density of a field ^{x) along 7 may be written in the form 

uVT,,(a;) = (T^Pi^*(v9 ® ^)) {x,x) (80) 

where the point split energy density operator is defined within a suitable neighbourhood 
W of 7 by 



a=0 



and {6^)0=0,1,2,3 is any smooth tetrad defined in a neighbourhood of 7 such that eg = u"" 
on 7. This operator may be used to define the renormalised energy density in the usual 
fashion. Given any real- valued / G C^{I), we may apply Theorem 13.11 in turn to the 
operators Qa = ■ V to obtain the absolute quantum weak energy inequality 



dr/^(r)(«VTr).(7(r)) > 



TC 



<,0 



+ [ dr/2(r)(g + MVC,6) 
Jr 



7(r) 



(821 



where : (r, r') 1-^ (7(r), 7(r')). 

For the purposes of comparison with existing QEI results, let us consider this bound 
for the massless Klein-Gordon field in Minkowski spacetime (M^, r]) for a worldline along 
the time axis. In this case, the bound simplifies because the full Hadamard series is given 
by the leading term; that is, A^{x,x') — l/(47r^cr+(x,a;')) is smooth and symmetric for 
any Hadamard state u. Of course, cr+ is globally defined in Minkowski space. Now 
consider the example above, applied to the case where 7 is an inertial curve. 



dr /2(r)KnXr)U7(r)) >-B'^- [ ^ 



(-^,0 (83) 
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(86) 
(87) 



where we have utihsed the fact that the Fourier transform of l/{t — iO^)^ is vr^^^(^)/3 
[3T] and changed variables to rj = ^ + C,. Hence, 

/ dtf(t)(«VTr).(7(t))>-T7^ / dr/l/(r/)|V (88) 
which is the same as the QEI for the massless field in Minkowski spacetime obtained in 

This example is of particular importance as on small length scales one expects the 
massive quantum field in a curved background to behave like its massless counterpart 
in fiat spacetime. We expect that the same should hold for the quantum inequalities, 
i.e., on small length scales the dominant contribution to the bound arises from the l/cr+ 
contribution to the Hadamard series. This will be investigated in a future work. 

3.2.3 QEIs for adiabatic states 

Finally, we show how our analysis of the Hadamard series using the Sobolev wave- 
front set allows us to establish QEIs for adiabatic states. We refer the reader to [3H] 
for a detailed study of adiabatic states and further references. Adiabatic states, like 
Hadamard states, are defined in terms of their singular structure: following [38jj^ a 
state uo on 2t(A^, (7) is an adiabatic state of order N if its associated two point function 
satisfies WF''{K^^) C for all s < + 3/2. From this definition, we see that any 
Hadamard state is adiabatic to all orders. In what follows the following lemma, taken 
from [38], is essential: 

Lemma 3.2. Let u be a Hadamard state and uo' be an adiabatic state of order N , with 
associated two point functions Ki^i respectively, on'Qi{J^,g). Then 

WF'{A^ - A^O = (89) 

for all s < N + 3/2. 

An immediate corollary is that ( l89i) also holds for all s < A^ + 3/2 if a; and u' are any 
two adiabatic states of order A^. By the Sobolev embedding theorem (Proposition 12.11 
(i)), differences of this type will be in C'^{Ai x A4) provided N > 9/2, thus permitting 
the construction of a normal ordered stress-energy tensor. Similarly, if u is adiabatic of 
order N > 9/2 and A; > 2, a difference of the form A^~Hk will be twice continuously dif- 
ferentiable on an ultra-regular domain, permitting the computation of the renormalised 
stress-energy tensor. It is straightforward to modify the proof of Theorem 13.11 to obtain 
the following. 

^In [35] the definition of adiabatic states was given, as for Hadamard states, only for quasi-free states, 
so the present usage is a shght extension. 
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Theorem 3.3. (a) Theorem \3.1\ continues to hold (with the same lower bound) under 
the weaker hypothesis that u is an adiabatic state of order N > 9/2. (b) Using the as- 
sumptions and notation of Theorem \3.1[ except that uj is assumed only to be an adiabatic 
state of order N > 9/2, there is a difference inequality 



dvol(x)/2(a;) Q®Q{A^- A^>) {x, x 



> -2 



h^\-^f^®\h^\~^f^dlQ®QK, 



'r+xK"-i (27r)" 

for any reference state uj' which is adiabatic of order N' > n + 11/2 



<,0, (90) 



Proof. We sketch the main points only. For (a), note that the hypotheses on N and k 
entail that we may choose s G (6, min{A^, A;} + 3/2). Part one of the proof of Theorem 13. II 
will continue to hold provided that A^ — Hj. G Hfoc- To see this, we introduce an arbitrary 
Hadamard state ujq and note that 

WF^A^ - Hk) C WF%A^ - A^J U WF%A^, - H^) = ^ (91) 

using Lemma r3. 21 together with the fact that s < N + 3/2, and Proposition 12.91 together 
with s < k + 3/2. Part two of the proof holds because A^^ is a bisolution to the Klein- 
Gordon equation, so all covectors in its wave-front set are null, from which it follows that 
the required pull-back exists. The third part is identical to the original argument. For 
(b), the hypotheses on N and N' permit us to choose s G (6, min{A^, N'} + 3/2), and we 
have WF^{A^ — A^/) = by the remark following Lemma Part one of the argument 
then goes through, as does the second part (with A^/ replacing Hk). The remaining 
issue is to check that the bound is finite. Introducing a reference Hadamard state ujq as 
before, we note that A^/ - A^^^, G H^^^i-^ ^ -^) some s' G (n + 7, A^' + 3/2), and 
hence A^/ — A^^ G C'^~^^{M. x M) (this is the reason for the constraint N' > n + 11/2). 
This is sufficient for Part three to apply to A^^i in place of Hj.. □ 



3.3 Covariance 

The QEIs obtained from theorem 13.11 are not covariant in full generality because they 
depend non-trivially on the coordinates used, and, in some cases, on a choice of tetrad 
near S. However, covariance may be rescued if we can restrict the freedom to choose 
coordinates and the tetrad in a covariant fashion so that the bound is independent of any 
residual choice. This strategy was successfully employed in [20] for worldline difference 
QEIs, and the same techniques would also apply to our worldline bounds; here we show 
how this may be accomplished for worldvolume averages (timelike submanifolds of other 
dimensions could be handled in an analogous fashion, but we will not do this here for 
brevity) . 

Consider the quantum null energy inequality studied in section 13.2.11 We will show 
that if the sampling tensor picks out a preferred smooth timelike curve 7 in a 
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covariant way and we employ a system of Fermi normal coordinates near 7, then residual 
choices in our construction cannot affect the bound. With these restrictions, our absolute 
QEI would be locally covariant in the sense of [20]; see also [23] for a more abstract 
discussion of these ideas in the formulation of locally covariant quantum field theory 
developed by Brunetti, Fredenhagen & Verch [1] in terms of category theory. 

The requirement that F"''' should pick out a unique timelike curve may be addressed 
in various ways. For example, if we restrict to sampling tensors for which there exists 
a (necessarily unique) pair of points x,x' E such that the support of the sampling 
tensor obeys 

supp F'^'' = J"(x) n J+(x') (92) 

and is contained within a convex normal neighbourhood then the unique timelike geodesic 
between x' and x may be used as our choice of 7. From now on we assume that the 
sampling tensor does indeed select a preferred timelike curve, and that 7 is given in a 
proper time parameterisation. 

As already mentioned, we will restrict our coordinate system to belong to the class of 
Fermi normal coordinates about 7. For completeness, we briefly summarise the salient 
features of Fermi-Walker transport and Fermi normal coordinates, mainly following 
chapters 1 §4 and 2 §10 of [46j- Recall that a vector field ^ defined on 7 is said to be 
Fermi-Walker transported along it if Dpw^ = 0, where 

dfwC = (7 ■ v)r - Qbcircy" - a^rje' (93) 

and a" = 7 ■ V7". Since a ■ 7 = and 7^ = 1 it is easy to see that Opwi = 0; hence, the 
velocity vector is preserved under Fermi- Walker transport. Moreover, it is possible to 
show that Fermi- Walker transport of two vectors along 7 preserves their inner-product. 
Therefore, a tetrad remains an orthogonal frame along 7 under Fermi- Walker transport. 
If 7 is a timelike geodesic, then Fermi-Walker and parallel transport coincide. 

The construction of Fermi normal coordinates near 7 proceeds as follows. Let y lie 
on 7 and construct an oriented and time-oriented orthonormal frame {e^}a=o,i,2,3 at y 
with Cq = 7" It/- Fermi- Walker transport yields a tetrad along the whole of 7. In a 
convex normal neighbourhood W of 7 each point x G W will be joined to 7 by a unique 
spacelike geodesic segment c which is orthogonal to 7 and which meets it at some 7(t). 
Assuming that c is parameterised by proper length, the Fermi normal coordinates of 
X are 

x° = t; x' = sc ■ ei\^(t) , (94) 

where s is the proper length of c. 

This construction has two important features. First, the metric takes the Minkowski 
form in these coordinates everywhere on 7. By continuity, this guarantees that the Fermi 
normal coordinates form a hyperbolic chart in a neighbourhood of 7. Second, the only 
freedom in the construction is the choice of origin on 7 (which amounts to the freedom 
to add a constant to x°) and the choice of the spatial tetrad vectors Cj {i = 1, 2, 3) at 
y, which are determined only up to a rotation. Owing to the angle-preserving nature of 
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Fermi-Walker transport, any two coordinate systems obtained by the construction are 
therefore globally related by = + X, x'* = S'^jX^ for constant scalar A and constant 
rotation matrix S G SO (3). 

If that the sampling tensor is supported within the neighbourhood of 7 in which the 
Fermi normal coordinates are hyperbolic, it is easy to see that the absolute QEI fl7^ is 
independent of the particular system of Fermi normal coordinates chosen. The key point 
is that the Jacobian determinant for a change of coordinates between two Fermi normal 
coordinate systems is identically unity by the remarks given above. (Note also that F"'' 
may be written uniquely as F"-^ = n"'n^ under the constraint that n"" is future-pointing 
and null where it is nonzero.) 

For more general QEls one also needs to construct a tetrad throughout the support 
of the sampling tensor. This may be done by taking the tetrad formed along 7 and 
propagating it by parallel transport along spacelike geodesies which meet 7 orthogonally. 

4 Conclusion 

We have given the first explicit absolute quantum energy inequalities for the massive 
minimally coupled Klein-Gordon field in arbitrary four- dimensional globally hyperbolic 
backgrounds, by refining the argument of [TU] to make use of the theory of the Sobolev 
wave-front set, and analysing microlocal properties of the components of the Hadamard 
series. The lower bounds are given in terms of partial sums of the Hadamard series, 
which are computed locally. Previously explicit absolute quantum energy inequalities 
were known only for the massless field in two dimensions [25] (a similar argument could 
be used to extend this to general positive energy unitary conformal field theories, based 
on [IE])- Although the bounds make use of coordinate systems, we have shown that by 
restricting the class of sampling tensors, there are circumstances in which the bound is 
CO variant. 

Absolute QEIs may also be found for higher spin fields. In the case of the Dirac field, 
which will be reported elsewhere [15] , one adapts the difference QEI obtained in [6] in a 
similar fashion to the way in which [TU] has been adapted here. Moreover, it is expected 
that one should be able to employ our method to prove an absolute QEI for the spin-1 
vector bosons, using the formulation of the Hadamard condition for the Maxwell and 
Proca fields given in p[^. 

One possibility which is opened up by our work is to obtain control over the size of 
spacetime region in which the absolute QEI bound can be approximated to a good degree 
by the QEIs obtained in Minkowski space for massless fields. This would potentially 
result in very simple bounds of wide applicability, and is the subject of ongoing work. 
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A Hadamard recursion relations 



In this appendix we briefly summarise the method for generating the coefficient functions, 
{vj}j=Q^,,,^k and {wj}j=o,...,fc; featuring in (1121) . These are obtained as the coefficients in 
the formal power series solution 

, 1 f A^x,x') ^ , a^{x,x') f a+{x,x') \ 

+ f:-^^M"-W^] ■ (95) 

to 

((V^ + fi^) ® l)Hk{x, x') = subject to wq = . (96) 

The series does not actually converge except in analytic spacetimes, which is why one 
makes use of the partial sums . The recursion relations for the Vj for the massive field 
in a curved background are: 

= f{\/'^ + fi^)A^+2Vvo-Va + Avo + voV^a (97) 

+ l)v,+i + {] + l)t;,+iVV (98) 

where j G {0} UN. In a regular domain X the system of differential equations uniquely 
determines the series of f/s. The Wj series is specified once the value of wq is fixed; we 
have adopted Wald's prescription that = [48j and with this boundary condition the 
recursion relations are: 

= 2Vwi ■ Va + WiVV + 2'S/vi ■ Vcr - Avi + t;i VV (99) 
= e{V^ + ii^)wk + 2{k + l)Vwk+i-Va 
-4:k{k + l)wk+i + ik + l)w;fc+iVV 

+2Vffe+i • V(T - 4(2A; + l)vk+i + f^+i VV (100) 

where A; G N. The system of equations (I97f98f99lll00p are known as the Hadamard 
recursion relations; these relations for the massless field may be found in PQ [7] where 
the dependency on a choice of length scale i is suppressed. 
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